Fully three dimensional breather solitons can be created using Feshbach resonance 
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We investigate the stability properties of breather solitons in a three-dimensional Bose-Einstein 
Condensate with Feshbach Resonance Management of the scattering length and confined only by a 
one dimensional optical lattice. We compare regions of stability in parameter space obtained from 
a fully 3D analysis with those from a quasi two-dimensional treatment. For moderate confinement 
we discover a new island of stability in the 3D case, not present in the quasi 2D treatment. Stable 
solutions from this region have nontrivial dynamics in the lattice direction, hence they describe fully 
3D breather solitons. We demonstrate these solutions in direct numerical simulations and outline 
a possible way of creating robust 3D solitons in experiments in a Bose Einstein Condensate in a 
one-dimensional lattice. We point other possible applications. 
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The creation of Bose-Einstein condensates (BEC) in 
vapors of alkali metals has opened an excellent opportu- 
nity to investigate nonlinear interactions of atomic mat- 
ter waves. An important challenge of practical interest 
is to develop methods to create and control matter-wave 
solitons. One dimensional dark 0, bright 0, and gap- 
mode solitons have already been observed. A promis- 
ing approach to obtain multidimensional solitons consists 
in varying the scattering length a of interatomic colli- 
sions. This can be achieved by means of sweeping an ex- 
ternal magnetic field through the a = point (the point, 
where the scattering length vanishes). This point occurs 
close to the Feshbach resonance 0. The application of 
an ac magnetic field may induce a periodic modulation of 
a, opening a way to "Feshbach-resonance management" 
(FRM) . A noteworthy FRM-induced effect is the pos- 
sibility of creating self-trapped oscillating BEC solitons 
(breathers) without an external trap in the 2D case. The 
underlying mechanism is fast modulations creating an ef- 
fective potential on a slower timescale. This potential can 
stabilize the soliton. The BEC model based on the Gross- 
Pitaevskii equation (CPE) with harmonic modulation of 
a was investigated in Refs. [6|, 121 13) • The conclusion was 
that FRM renders it possible to stabilize 2D breather 
solitons even without the use of an external trap. Ac- 
cording to these references, 3D breathers require at least 
a tight, one dimensional harmonic trap jfj, practically 
reducing the problem to 2D. Later on we will call this 
approach a quasi two dimensional (Q2D) treatment. It 
will be defined more precisely below. 

Recently flOj we demonstrated that quasi 2D solitons 
can be stabilized by a combination of FRM and a strong 
ID optical lattice (ID OL), instead of a ID harmonic trap 
jH EL @ ■ By a "strong lattice" we mean one in which the 
atoms in neighboring cells cannot interact. This issue has 
practical relevance, as a ID OL can easily be created, illu- 



minating the BEC by a pair of counterpropaga ting laser 
beams that form a periodic interference pattern . The 
lattice will be weak or strong depending on the intensity 
of the laser light. In fact, it is easier to realize a tight 
confinement configuration in an optical lattice than in 
an harmonic trap. Hence this environment may be more 
friendly for creating quasi 2D solitons. 

In this Letter, we demonstrate that the combined OL- 
FRM stabilization of 3D solitons is possible even in a 
weak lattice, when atoms confined in different cells in- 
teract. By analyzing the stability charts in configuration 
space we discover two distinct regions where stable solu- 
tions exist. The first of these regions has its counterpart 
in the Q2D treatment. The other region appears when 
the frequency of modulation exceeds a critical value de- 
pendent on the confining potential. It is not present in 
the Q2D treatment; hence it corresponds to fully 3D soli- 
tons. In the limit of tight confinement the latter region 
moves to extremely high frequencies and the Q2D stabil- 
ity chart is recovered. 

We describe our system by the GPE in physical units, 
including a time-dependent (FRM-controlled) scattering 
length a(r) and an external potential U(r, t) 
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Initially the BEC is in the ground state of a radial (2D) 
parabolic trap with frequency <Dj_, supplemented, in the 
longitudinal direction, by "end caps" induced by trans- 
verse light sheets. The configuration is much like the one 
used to create soliton trains in a Li 7 condensate 0- A 
ID lattice potential in the axial direction is adiabatically 
turned on from i = to i = if , see Fig.^ Thus, the full 
potential is 

U(r, t) = e(r) [1 - cos(2ttz/A)] + 
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+ f(r) (m/2)^ig 2 + U (z) 
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where A is the lattice spacing, g is the radial variable 
in the plane transverse to z, and the axial "end-cap" 
potential, Uq(z) is approximated by a sufficiently deep 
one dimensional rectangular potential well. The width 
of the well determines the number of peaks in the fi- 
nally established structure. The /(r) is a switching-off 
function (see Fig. ^) . We introduce dimensionless vari- 
ables x = (7r/A)r, t — tuj, ip = ^ ^ X 3 / (Nn 3 ), where 
u) = Tr 2 h/ (mX 2 ) and N is the number of atoms. Wc 
also define g = 4n 2 Na/X, e — e/(uK), uj± = Oj±/uj, and 
Uq = Uq/(u>K) to obtain 
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and the potential in the form 
t/(x, t) = e(t) [1 - cos(2z)]+/(<) 
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The nonlinear interaction coupling is described by 
9(t) =9o(t) +gi(t) sin(Ot), 



(3) 



(4) 



(5) 



and the dimensionless modulation frequency by f2 = 
Cl/u. Initially gi (0) = and 5(0) = g (0) > 0. At some 
moment t%, we begin to linearly decrease go(t). It van- 
ishes at time £2, and remains zero up to £3, when we start 
to gradually switch on the rapid FRM modulation of a. 
In the interval [£3, ti] , go(t) decreases linearly from zero to 
a negative got and the amplitude of the modulation g\ (t) 
increases from zero to gu , see Fig. ^ Simultaneously, 
both the radial confinement and end-caps are gradually 
switched off. At times t > t 4 , g(t) oscillates with a con- 
stant amplitude gu around a negative average value got. 
Consequently, a soliton so created, if any is supported 
by the combination of the ID lattice and FRM. 

Numerical experiments following the path outlined in 
Fig. 2] indicate that it is possible to create stable solitons 
[lfj (see inset to Fig|2J). Before showing the results, we 
first resort to the variational approximation (VA) in order 
to predict conditions on the modulation frequency and 
the size of the negative average nonlinear coefficient got , 
necessary to support 3D solitons. 

The VA can be applied to the description of BEC dy- 
namics under diverse circumstances H E H E El 
Equation Q is derived from the Lagrangian density 

c = i(r t Tp - - - iv^i 2 - gWM 4 - 2u\^\ 2 (6) 

We use VA for t > t% and choose a complex Gaussian 
ansatz for the solution for one lattice cell. The amplitude 
is A(t), radial and axial widths are W(t) and V(t) re- 
spectively, and b(t) and (3{t) are the corresponding chirps 

ip(r, t) = J 4 e [-e 2 (i/2W 2 + l b)-^(i/2y 2 + l/ 3)+ 4 0] _ ^ 
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FIG. 1: (Color online). The time dependence of the nonlinear 
coefficient, g, switching-off function, f(t), and optical-lattice 
strength, e, in the numerical experiment which leads to the es- 
tablishment of stable 3D breathing solitons supported by the 
combination of the quasi-lD lattice and Feshbach-resonance 
management (FRM). The shaded area indicates rapid oscil- 
lations of g(t), which account for the FRM. 



The reduced Lagrangian can be found upon substitut- 
ing J7J into ijBJ) and integrating over space. By varying 
this reduced Lagrangian with respect to 4> we obtain the 
constant E = A 2 W 2 V = tt" 3 / 2 J cell \*P\ 2 dr = ir~ 3 / 2 /n, 
where the integral extends over one cell of the lattice and 
n is the number of occupied lattice cells. Notice that the 
total number of atoms is included in the definition of the 
nonlinear coupling g(t) and the total wavefunction is nor- 
malized to unity. When the other four variational equa- 
tions are derived, we can deduce two dynamical equations 
for the widths: 
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(9) 



These equations describe the dynamics of a single peak, 
and with one term slightly altered can be applied to the 
problem of a BEC confined in a ID harmonic trap 0, 0, 

a 

In the corresponding Q2D treatment we drop the z 
dimension in Eq. @. It is assumed that in this direc- 
tion the profile of the wavefunction is fixed and repro- 
duces the ground state tpo of the single lattice cell (Wan- 
nier function) or harmonic potential, as in Refs. 
The reduced potential in 2D GP will take the form 
U(g,t) = f{t)\l/2)uj 2 _g 2 . In VA we take V = V Q from 
4erfV 4 exp (— V 2 ) = 1, and only solve equation JSJ. In 
numerical simulations we rescale the nonlinear coupling 
coefficient <?2D = g X {J \i>o\ 2, 4>odz)/(f ipodz). 

We simulated both the full GPE, Eq. using an 
axisymmetric code (for 3D), a Cartesian code (for 2D), 
and the variational equations for comparison. Numerical 
simulations followed the path outlined in Fig.^ The pa- 
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FIG. 2: Evolution of the amplitude of the central peak in 
the pattern comprising five cells of the optical lattice (shown 
in the inset). The normalized parameters are gof = —30, 
gi{ = 4p f, £f = 20.5, Q = 22, lo± = 0.3, ti = 30, t 2 = 100, 
ts — 120, and tt = 130. The snapshots in the inset is taken 
at t = 6000. The unit of time is m\ 2 /{-R 2 h). 




rameters used in the simulations would correspond, for 
85 Rb atoms and an OL period of A = 3 /im, to an initial 
radial-confinement frequency of uj± = 2ir x 39 Hz, a FRM 
frequency of = 2ir x 2.87 kHz. A good candidate is the 
Feshbach resonance at 155 G for 85 Rb. Since in this case 
the a = point is far from the resonance, both two and 
three body losses are negligible 14]. A lattice depth of 
£f = 10.25 recoil energies, and an effective nonlinear co- 
efficient of {N/n)a — ± 10 -5 m are required. Here N/n is 
the number of atoms per lattice cell, with a total number 
of atoms in the range of 10 4 — 10 6 . The respective val- 
ues of the normalized parameters are given in the figure 
captions. 

Examples of the numerical results, which are generic, 
are displayed in Fig. [21 which shows the evolution of 
the central-peak's amplitude versus time in dimcnsionlcss 
units, defined by Eq. ©. After an initial transient, a sta- 
ble structure is established, featuring breathing without 
any systematic decay. The figure reveals the influence 
of neighboring solitons on the amplitude of the central 
peak. The thin curve corresponds to the evolution of 
the full multipeak structure, shown in the inset for a 
fixed moment of time. To obtain the thick curve we re- 
peated the above calculations up to the time t = 7500 and 
then removed all but the central peak. In this case the 
amplitude gradually decreased. The interaction between 
neighboring solitons can be explained by the difference 
between oscillation periods from cell to cell due to small 
deviations in the numbers of atoms. It is essentially the 
same mechanism as in a Josephson junction. When the 
lattice potential is weak, interaction between neighboring 
solitons is possible. 

In Figs. |21 and [I] we have collected results of a sys- 
tematic scan of parameter space based on GPE simu- 



FIG. 3: (Color online). Stability regions for the 3D solitons 
in the (|ffof ]) ^) plane, as predicted by the variational approx- 
imation (shaded area), and found from direct simulations of 
the Gross-Pitaevskii equation (circles) . Other parameters are 
as in Fig. |5] except Sf = 50. Note the similarity of the lower 
region to that of Fig. H This region corresponds to Q2D 
solitons. On the other hand, the upper region contains fully 
3D solitons. The borders of the VA stability regions were 
found analytically and will be given in a fuller version (l5|. 
The lowest excitation frequency of the confining potential is 
fio = 26.76. 




FIG. 4: (Color online). Same as Fig. 03 but for Q2D treat- 
ment. 
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lations and compared them with the predictions of the 
VA (a similar analysis can be performed if we replace ID 
OL with a ID harmonic trap - the conclusion does not 
depend on the form of confinement in the z direction). 
By stability we mean shape preservation during one run. 
When this was not the case, after a short period of time 
we observed clear collapse or spreadout in transverse di- 
rections. 

The agreement between VA and direct simulations is 
very good. As seen from Fig. EI m °ur fully 3D treatment 
we found two islands of stability. Note the similarity of 
the lower region to that of Fig. which portrays the re- 
sults of the Q2D treatment. This region corresponds to 
Q2D solitons. On the other hand, the upper region con- 
tains fully 3D solitons. It appears when the frequency 
of modulation exceeds the lowest excitation frequency of 
the confining potential, Slo = 26.76. If the strength of the 
lattice £f is increased, this region moves towards higher 
frequencies, the Q2D region expands upwards, and be- 
comes more and more like in Fig.0J This will be demon- 
strated in a fuller version of this work 0] . 

As we saw, new region of stability appears in the 3D 
treatment as compared to the Q2D treatment. We are 
more used to the effect of adding a new dimension simply 
shrinking or abolishing the basin of stability of solitons 
or waves. For example, water waves are unstable with 
respect to perturbations along their direction of propa- 
gation only when the depth exceeds a critical value |16| . 
When, however, two dimensional perturbations are al- 
lowed, there will always be an unstable angle regardless 
of the depth ^3] • Another example is that of ID solitons 
of the Nonlinear Schrodinger equation (NLS) with con- 
stant coefficients. They are stable in ID, but unstable in 
2D or 3D. This is also true for some NLS waves |l8| . 

In the problem treated in this Letter this is not the 
case. For much of the quasi two dimensional stability 
chart adding a degree of freedom stabilizes the soliton 
solution. The key to this dichotomy seems to be the 
presence of a periodic modulation, absent in the above 
mentioned classical examples. This can be illustrated 
by a simple case involving oscillators. Take as the one 
dimensional version a forced oscillator problem: 

x + ujqX — y cos(woi). (10) 

If y is fixed, the solution has a secular component x = 
yt/(2uJo) sin(wo^) + F{t), where F(t) is a periodic func- 
tion, and so the amplitude will grow as t. If however we 
allow a second degree of freedom, such that y also oscil- 
lates (for instance y + e 2 y = 0) the solution stabilizes, 
unless e = ±2uiq. In general this can be the case when 
there are periodic modulations. This fact, obvious in os- 
cillator theory, is perhaps less well known in the soliton 
context. 

The main result of this letter is the possibility of cre- 



ating fully 3D breather solitons in a BEC confined by a 
ID optical lattice potential, corresponding to the upper 
region in Fig.|3 The stable patterns may feature a multi- 
cell structure, which in the case studied here forms a set 
of weakly interacting fundamental solitons. The scheme 
proposed in this Letter is based on a combination of FRM 
and a ID optical lattice, and could be implemented in an 
experiment, as outlined here. This would open the way 
to the creation of robust 3D solitons (breathers) in BECs. 
A similar idea can apply in the field of nonlinear optics. 
The Feshbach resonance could be replaced by a nonlinear 
periodic structure, for example ferromagnetic domains as 
used in contemporary quasi phase matching . or by 
optically induced photorefractive lattices |2(| . 
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